POLYNOMIAL DIFFERENCES IN THE PRIMES 



NEIL LYALL ALEX RICE 



Abstract. We establish, utilizing the Hardy-Littlewood Circle Method, an asymptotic formula for the 
number of pairs of primes whose differences lie in the image of a fixed polynomial. We also include a 
generalization of this result where differences are replaced with any integer linear combination of two primes. 



1. Introduction 

Given a natural number N, how many pairs of primes less than or equal to ./V differ by a perfect square? 
More generally, for an arbitrary polynomial / £ Z[x], we define 
(1) r f (N) = #{{pup 2 ) e P% : Pi -P2 e /(N)}, 

where V denotes the primes and Vn — V (1 {1, . . . , N}. In this article we will first provide a heuristic, then 
argue rigorously with the Hardy-Littlewood Circle Method, for the following result. 

Theorem 1. If f(x) — CkX k + ■ • • + C\X + cq € Z[x], c% > 0, and k > 1, then 

^Jp\ (P - I) 2 J cl /k k + 1 log 2 N \ log N 

where Zf(p) denotes the number of roots of f modulo p, and the implied constant depends only on f . 

We note that due to the sufficient decay of (z/(p) — l)/(p— l) 2 , the product in the leading term is always 
finite, and is only if one of the individual terms is 0. This only occurs if z/(2) = 0, yielding a meaningful 
result as long as / has a root modulo 2. Also, if f(x) = cx k , then the product collapses to Yi P \ c ^T' including 
its complete disappearance if f(x) = x k . 

Remark. The condition that / has a root modulo 2 is one piece of the more stringent condition that / is 
an "intersective polynomial" , i.e. / has a root modulo n for all n G N. This additional assumption on / 
is necessary and sufficient to conclude the same order of growth for the count analogous to rt(N) for any 
subset of the primes of positive relative upper density, shown in a recent result of Thai Hoang Le [4] . 

While there is an extensive literature devoted to questions of this type going back some 60 or 70 years, 
namely Goldbach type problems for polynomials (see for example [3], QZI], [Z] and [T]), it appears that the 
precise result stated in Theorem [1] above has not been considered previously. 

2. Heuristics 

We begin by addressing the original motivating question, temporarily letting r(N) = r x 2(N). If we 
consider the discrete derivative r'(N) = r(N) — r(N — I), we see that r'(N) = unless N g V, which by the 
Prime Number Theorem occurs with probability 1/logiV. If N does happen to be prime, then it induces 
I'Pzv-il ~ N/ log N new positive differences of primes, which are each perfect squares with probability 1/V~N. 

Therefore, we may well expect r'(N) to grow on average like y/N/ log 2 N, which in turns leads, via a 
simple summation by parts argument, to the prediction that 

N ? A/3/2 

(3) r(JV)=yV(fc) = - 
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This heuristic immediately generalizes to predict 

(4) r x k(N) = ^ h error, 

W x K ' k+l \ og 2 N 

but it fails to account for potential congruence biases. 

In the spirit of generalization, we fix f(x) = c^x k + ■ ■ • + C\X + c G Z[x]. Note that rj(N) — r_f(N), 
so we can assume without loss of generality that c,t > 0. We must consider how rt(N) should compare to 
the prediction ^ above. One difference is that the leading coefficient Ck causes the number of points in the 
image of / that are at most N to be reduced asymptotically by a factor of c k . We must also consider local 
congruence biases of both the primes and the polynomial. For example, if we took f(x) = 2x k + 1, then the 
second prime in each pair would have to be 2, rt(N) would clearly then be bounded by \Vn\ ~ N/\ogN, 
and we couldn't possibly obtain the same order of growth as predicted in (|4]). The reason for this collapse 
is of course that almost all primes are odd, so almost no differences of primes are odd. More generally, we 
can use knowledge about the distribution of primes to gain knowledge about the distribution of differences 
of primes, and use these facts to investigate which polynomials these differences favor and avoid. 

In order to exhaust all possible congruence biases, we must consider the distribution of differences of 
primes and the image of our polynomial modulo n for each n £ N. However, if two moduli n and m are 
coprime, then congruence modulo n and congruence modulo m are independent events (by the Chinese 
Remainder Theorem). Therefore, we can determine all bias by investigating modulo arbitrarily large powers 
of each prime. Toward this end, we consider an arbitrary prime p, and we proceed probabilistically. We note 
that if two sets A and B are uniformly distributed modulo p n , and we randomly select a € A and b £ B, 
then the probability that a = b modulo p n , which we denote by P(a = b mod p n ), is l/p n . We compare the 
analogous probability for our specific sets with this expectation by defining a "bias factor" 

(5) b f (p n ) = P(Pi-P2 = /W™dp") 
and further noting that 

p«-i 

b f (p) = p n J2 P (Pi ~P2 = a mod p n )P{f(d) = a mod p n ) 

a=0 
P n -1 

= P(pi — P2 = a mod p™)(# solutions to f(x) — a in Z/p n Z). 



Once we determine these biases, we can make a prediction of the form 

(6) r f (N) = [] limsupfc^) — + error. 

p 6 -p c k ' K + 1 log N 

This formulation turns out to be unnecessarily frightening. The Prime Number Theorem for Arithmetic 
Progressions due to Siegel and Walfisz states that for any modulus m, the primes are evenly distributed 
among the congruence classes coprime to m. Therefore, the only biases of the primes, and hence the only 
biases of differences of primes, are related the coprimality. However, we know that (a,p) — 1 if and only if 
(a,p n ) = 1 for all n, which means that P(pi — P2 = a mod p n ) = P{pi — P2 = b mod p n ) if a = b mod p. 
One can easily show from the definition above that this implies bf(p n ) — bf(p) for all n, greatly simplifying 
prediction (j6|) to 



(7) r f (N)=Y[b f (p)4 / , ] ±- i 



N (k+i)/k 
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error. 



Now, for each fixed prime p, it follows from the Siegel- Walfisz Theorem that primes are congruent to 
1, 2, ... ,p — 1 modulo p each with probability l/(p— 1) and consequently that the difference of two randomly 
selected primes is congruent to modulo p with probability l/(p — 1), and congruent to 1,2, ...,p— 1 
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modulo p each with probability (p — 2)/(p— l) 2 . From this observation it follows that the bias factor bf(p) is 
completely determined by the number of roots of / modulo p, a quantity we shall denote by Zf(p). In fact, 

( 8 ) b f(p) = ^J z fiP) + (P ~ */(?)) = 1 + ^3T)F- 

This completes the heuristic for the result in Theorem [T] and indicates that the only "fatal" obstruction 
toward the expected order of growth for r/(N) is the modulo 2 consideration noted above. 

3. Rigorous Treatment via the Hardy-Littlewood Circle Method 

In this section we give a proof of Theorem [T] using the circle method developed by Hardy and Littlewood. 
As is standard, we begin by weighting the characteristic function of the primes with a logarithm to obtain a 
more uniform distribution. This yields a weighted count intimately related to r/(iV), which we define below 
using a truncated von Mangoldt function. 



\ogp if n — p k < N,p eV,k G N 
else 



M 



1 h { 
(11) R f (N) = &(f) —N^+W + O ^ 



Definition 1. We define Ajv : Z — > [0, oo) by 

(9) A N (n) 
and for f € Z[x] as in Theorem 1, we define 

(10) i? / (A)=^^A JV (n)A JV (n-/(d)), 

d=l nSZ 

where M = (N/c k ) 1/k . 

The main result of this section, and indeed the whole article, is the following. 
Theorem 2. For f G Z[x] as in Theorem 1 and any A > 0, we have 

c k 

where the implied constant depends only on f and A, and 

(12) 6(/)=E^p E E^ /(r)a/9 ' 

q=l 0<a<q r=Q 

(a,q) = l 

where fi is the Mobius function, (j) is the Euler totient function. 

It is easy to see that Theorem [1] will follow from Theorem [2] and the following two lemmas. 
Lemma 1. If f G Z[x] is as in Theorem\]\ and Zf{p) denotes the number of roots of f modulo p, then 

< 13 > £ i:w = n(i + ^)- 

g=l ^r\Ht < a<9 r=0 pgP v \f 1 / 

(a,q)=l 

Lemma 2. If f € Z[x] is as in Theorem^ then 

(u) rf{N) = ^m +0 (^yyio S N 

' n ' log N V log N 

The proof of both Lemma [T] and Lemma [2] are standard exercises, for completeness however we include 
their proofs in Section |4] below. In the remainder of this section we present the proof of Theorem [2] 
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3.1. Proof of Theorem [2j We fix / as in Theorem [T] and A > 0. 

In order to apply the circle method, we need to reformulate the weighted count Rf(N) on the transform 
side as an integral over the circle. Specifically, it follows from the usual orthogonality relation 



(15) / e 2ma da 

that 



1 if n = 
if n G Z \ {0} 



p 2tt ina 



1 

(16) R f (N) = I \A^(a)\ 2 S M (a)da, 
where 

(17) A^(a) = Y, A N{n) 

nez 

and 

Af 

(18) S M {a)=Y. e2 * lf[d)a - 

d=l 

As usual, we wish to partition the circle into two pieces, one the collection of points near rationals with 
small denominator, from which the primary contribution to the integral (|16[) will stem, and the other simply 
the complement. In order to appropriately make this partition official, we need a parameter yielded from 
the following estimate on high moments of the Weyl sum Sm ■ 

Lemma 3. There exists so(fc) = 0(fc 2 logfc) such that if s > so(k), then 

(19) [ \S M (a)\ s da = 0(M s - k ). 

Jo 

For a proof of this estimate see Proposition 3.3 in j5]. We remark that it follows from recent work of 
Trevor Wooley [5], that one can in fact take So(fc) = 2k(k + 1). 

Definition 2. Fixing s = max{so(k), 4fc} we let B = s(A + 1) + 8 and define 9Jt, the major arcs, by 

log s N 

(20) 971= |J |J M a/q , 



9 =1 0<a<q 
0,S)=1 



where 

(21) M a/ , = L[fl,l] 

and m, the minor arcs, by m = [0, 1] \ 971. 



a 

a 
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N 



3.1.1. Estimates on the minor arcs. We first focus our attention on the minor arcs, in an attempt to absorb 
their contribution to the integral into the error term in Theorem [2] To this end, we invoke the minor 
arc estimate on Ajv that was Vinogradov's main achievement in solving the ternary Goldbach problem 
unconditionally. 

Lemma 4 (Vinogradov, Vaughan). If a £ [0, 1] , \a — a/q\ < 1/q 2 and (a,q) = 1, then 
(22) 7W(a) = O (log 4 N (N/q 1 ' 2 + N^ 5 + N^q 1 / 2 ^ ) . 

The argument for this estimate was greatly simplified by Vaughan, and an exposition of it can be found 
in [8], for example. It is now a relatively straightforward matter to see that Lemmas [3] and |4] combine to 
give the following desired estimate for our integral over the minor arcs. 
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Corollary 1 (Minor arc estimate). 

' M(k+l)/k 

(23) / \A N (a)\ 2 S M (a)da = 0' 



\og A N 



Proof. First we fix an arbitrary a € m. By the pigeonhole principle, there must exist 1 < q < N / log 5 N 
such that 



a 

q 



log" N 1 
< — < — . 

qN ~ q 2 



However, by the definition of m, it must be the case that q > log" 8 N. Combining these bounds on q with 
Lemma |4] gives 

(24) A N (a) =0' A 



Now, we invoke Holder's Inequality in order to utilize our control on the higher moments of Sm ■ Indeed, 
applying Holder's inequality, Lemma [3] and (|24p we obtain 

\A N (a)\ 2 S M (a) da 



< (J \A N (a)\ 2s ^ s -V da\ (J \S M (a)\ s dc^j 

^((w^) 2/s (r^ (a)|2da ) (s 1,7 m " 



o 



/ 7v(fc+i)/fc 
V log A N 



In the last line above we have used the fact that M = 0{N 1 / k ) and B — s{A + 1) + 8, together with the 
deeper fact that 



l 

2 



\A N (a)\ 2 da = O(NlogN), 

which is a standard consequence of the Prime Number Theorem (and Plancherel's Identity). □ 

3.1.2. Estimates on the major arcs. For the major arcs, we invoke the most classical of estimates. 
Lemma 5 (Major arc estimates). If a = a/q + /3 with 1 < q < log 5 N , (a, q) — 1 and |/3| < log 5 N/N, then 
(25) S M (a) = q- X S a/q I M ((3) + 0(\og 2B N) 



(/// 



d 



(26) A N (a) = ^vnW) + 0(7Ve- c ^ W ) 

for some c — c(B) > 0, where \x is the Mobius function, (f> is the Euler totient function, 

Sa/ q = Y, e27Tinr)a/q > M£>=/ e^^dx, and v N {[3) = ^ e" 

r=0 ^° n=l 



27rm/3 



Estimate (|25|) can be found in any discussion of Waring's Problem (Lemma 4.2 of [2] for example), while 
estimate ([26")) follows from the Siegel-Walfisz Theorem on primes in arithmetic progressions, and is included 
in any discussion of Vinogradov's Three Primes Theorem (see for example Lemma 3.1 in [5]). We further 
note that 

(27) v N (fi)= / e- 2mx0 dx + O{\og B N), 

Jo 

thus we maintain the same asymptotic by replacing the sum with the integral. 

We are now ready to combine all of our estimates and establish Theorem [2l 
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3.1.3. Proof of Theorem^ From (j!6[) and Corollary [TJ and noting that for large N the individual major 
arcs are pairwise disjoint, we have 

l °s B N r /\T(k+l)/k\ 

(28) Rf(N)= J2 E / \^N(a)\ 2 S M (a)da + 0[ A — . 

(o,?)=l 

From (|2"5]l. (USD, and (|27]l. we then obtain 



log" AT 



(29) %(A0= £ 



=1 ^(<?) 2 ^ <q 



(aT,)=l l fl<(log B iVVf 







/ 


/ e " 27ri ^ dy 




J 


; B iV)/JV 





I M {P)dp+0 



N (k+\)/k 
\og A N 



(*) 



(**) 



The reader can now easily see that Theorem [5] will be an immediate consequence of the following estimates 
for (*) and (★★) respectively: 



log" N 



(30) 

and 
(31) 



U ot< q 

(a,q)=l 



i r 



-27riyj3 



dy 



|/3|<(log B N)/N 



(o,g)=l 



I M (l3)dl3 



k+1 



MN + 



los A N 



MN 
log 8 N 



In order to establish (j3"0")) we need an estimate on the magnitude of the Gauss sum S a / q that beats the 
trivial bound of q. Theorem 7.1 in [5] tells us that 

(32) S a/q = O e (g( 1 - 1 / fc )+ e ) for all e > 0, 

so in particular, we know that S a/q = Oiq 1 - 1 /^), and since <j)(q) > Cq 1 ^ (trivially), we can deduce from 
this that 



E 



M<z) 2 



(o,?)=l 



E S «/9 - I E a l + l/4fc 
<a<9 y 9 =log B AT ^ / 



l 0K B/4fc iV 



Noting that B/Ak > A, we see that (l3"0"|) immediately follows. 
In order to establish (pTTj) we initially note that 



a/ 



M0) - / e 2 " Cfea: ~ p dx = 0\ x k - x pdx = Oi\og B N) 



(33) 



and that substituting this into the left hand side of pip gives 







/ 


/ e - 2nyl3 dy 




Jo 


|/3|<(log B N)/N 





hi(P)d(3 = 

|/3|<(log B N)/N 



N 



dy 



e 2wic * xk Pdx\ d/3+0(N log 25 N). 



We note further that after three changes of variables (namely x := x/M, (3 := N/3, y := y/N), the integral 
on the right hand side of the above identity can be seen to satisfy 
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(34) 



|/3|<log B N 



-2iriyf3 



dy 



M 



dx\ dp = MN / llj^^)! 2 / e 2 ™ k Pdx) df3, 



|/3|<log B N 



where l[o,i] (/?) = Jo e 27 ™^ dy is the usual Euclidean Fourier transform of the characteristic function of the 
unit interval. Utilizing the standard (and easily verified) fact that |l[o,i](/3)| 2 = 0(1 + |/3| 2 ) -1 , it then follows 



that 

(35) J |1^J(/3)| 5 

|/3|>log B N 

and hence that 



C e 2mx ^dx) d(3 = o( —tr- 
io ) \\og B N 



(36) J IW/?)| 2 Qf^-^dx) d(3 = Jj^m 2 Qf* e 2 "**^ df3 + O (r^) • 

|/3|<log B JV 

Thus, in order to establish (j3"Tj) it finally suffices to show 

(37) jf |f^(/3)| 2 Qf 1 e 2 ™*^*) ^ = fcTT- 

Changing variables again (a: := £ fc ), letting gr(ar) = x^~ k ^ k /k l[o,i](a;), and applying Parseval's Identity, 
it follows that 



l[o,i]0*0 / ff(y)l[o,i](s-!/)<^da; 



1 r 2/ (i_fe)/fc 

/ — E — rf y dx 

Jo fc 
_ k + 1 

as required. □ 

4. Proof of Lemmas □] and [2] 
4.1. Proof of Lemma [TJ First we write 

(38) ESE £e-/0>/^ $>(«), 



g=l 0<a< 9 r-0 9=1 

(o,g)=l 



where 



/ \2 i- 1 
H'KQ) ST^ G „„j O _ e 2Ttif{r)a/q 



(39) F(g) = — T-rj £ 5 »/9 and 5 <V<? = 

(o,?)=l 

Since q, n(q), and </>(g) are all multiplicative functions, and it is a standard exercise to verify that 

(40) ]T s a/q 

0<a<q 
(a,q) = l 

is also multiplicative (see for example Lemma 5.1 of [2] for a proof), it follows that F(q) is a multiplicative 
function. Furthermore, during the proof of Theorem [2j we used estimate ff32]) to establish that F(q) = 
0(q~ 1 ~ 1 ' 4k ), so in particular we know that X^i m an absolutely convergent sum. 
These two facts combine to yield the usual Euler product formula, namely 

oo oo 

(41) e F( q ) = n e w = n o- + f(p)), 

9=1 pGV n=0 p£V 
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with the great simplification occurring due the presence of the Mobius function, which forces F(p n ) = for 
every prime p, provided n > 2. 

Noting that 



(42 ) ^) = ^^EE e27ri/(r)a/p ! 

' a=l r=o 

we see that it now suffices to simply verify that 

p— l p— l 

(43) Jj;^/W«/p=p( 2/ (p)-l) 

r=0 a=l 

for all primes p. This is straightforward, since if f(r) =0 modulo p, then the exponential term is identically 
1, and 

p-i 

e 27ri/(r)a/p =p _l 

a=l 

While if f(r) ^ modulo p, then the inner sum goes over all p-th roots of unity except 1, hence 

p-i 

e 2*if(T)a/p _ _j 

o=l 

Therefore, since z/(f>) denotes the number of roots of / mod p, it follows that 

p— l p— l 

(44) ^ ^ e^/O/P = ^(pXp - 1) - (p - z,(p)) = p( 2/ (p) - 1) 

r=Q a=l 

as required. □ 
4.2. Proof of Lemma [2j First, we let 

(45) M' = min{n e N | /(d) > N for all d > n}, 
so 

M ' / N 

( 46 ) EE 1 ^W 1 P»( ,l -^) = #fe d ) ePffXN : P-/(d)6P w }=r/(iV) + 0' 

d=l n£l, 



log iV / ' 



where the error term accounts for the finitely many instances of "double-counting" in the polynomial. This 
convenience is the primary reason for restricting to the image of N as opposed to Z, making all nonconstant 
polynomials "eventually injective". We also note that M' = M + O e (N e ) for every e > 0, and in particular, 



M 

l/4fe 



(47) r f (N) = E E Wn)Wn - /(d)) + O (N 1+1 



d=l raSZ 

We now make the usual observation that the contribution to the von Mangoldt function from proper prime 
powers is negligible. Namely, 

M 

(48) R f (N) = 53 53 log(n) log(n - f(d))l VN (n)l VN {n - /(d)) + O (nW™ log 2 n) . 

By (HI]), f4"7| . and ([IS]) , noting the trivial upper bound 

M A/ 

(49) 5353 l VN (n)l VN (n- /(d)) > — y- 53 53 log(n) log(n - /(d))l Pjv (n)l Pjv (n - /(d)), 



and changing variables (n := n — f(d)), it now suffices to show 

M M 

£ £ ^« ip)l VN (n + /(d)) < — J- £ £ l0 S( n ) lo §( n + (n)lv N (n + /(d)) 

/ jy(fc+i)/fcl oglog j y- 
V log 3 TV 
Toward this end, we see that for any <5 > 0, 

M M N M N 1 - 6 

^^lp N (n)l PN (« + /(d))=^ Y lvAn)lvAn + f(d))+Y £ lv N {n)lv N {n + }{d)) 

d=lneZ d=l n =Af 1 -' 5 d=l n=l 

M 

< £ £ log(n) log(n + /(d)) Wn)!^ (n + /(d)) 



+ 
Setting 



(l-^log 2 ^; 



\N S log TV 



5= 2 log log iV 



log TV ' 

and noting from Theorem [5] and (|48l) that 

M 

( 51 ) £ £ lQ g(") lo §( n + / W) 1 ^ (" + /(<9) = (iV (fc+1)/fc N 

<i=l riGZ 

the lemma follows. □ 

5. Further Generalization 

A natural question: How would our results change if, instead of differences of primes, we counted sums of 
primes that lie in the image of a fixed polynomial? Even more generally, for / £ Z[x], TV G N, and arbitrary 
integers ai,a 2 ^ 0, we define 

(52) r /l01 , 0B (JV) = #{(p 1 ,p 2 )e^ : a lPl + a 2 p 2 G f(N)}. 

Again we note that r/ l0l) a 2 (TV) = r_y._ ai ._ a2 (TV), so we can assume that / has positive leading coefficient. 
Under this assumption, if ai and a 2 were both negative, one can see that rt iaiia2 would be uniformly bounded 
in TV, so we assume without loss of generality that a% > 0. 

Careful adaptation of our initial heuristic with the discrete derivative yields an additional factor of 

f(ai+a2) 1/fe -a 1 /fc (a x + <x 2 ) 1/fc - a 2 /k 



(53) 



if a 2 > 

a 2 a\ 
ni \ I (ai + a 2 ) 1/k - a\ /k (ai+a 2 ) 1/fc 

gfc(ai,a 2 ) = S y — + — if0>a 2 >-ai 

a 2 ai 

l/fc 

if a 2 < — ai 



n 2 

Note for example that Cfc (1,— 1) = f for all k, as it should to agree with our previous discussion, and 
C fc (l,l) = 2(2^ -1). 

To adapt our heuristic for local congruence biases, it is useful to partition the primes based on how many 
of the coefficients a\,a 2 they divide. Namely, we define 

Vo = {peV :p\ ai a 2 } 

(54) Vi = {p G V : v I aia 2 , (p f ai or p f a 2 )} 

^2 = {p G T 5 : p | ai, p | a 2 }. 
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Hopefully context will prevent any confusion between this notation and the notation Vn = V P\ {1, . . . , A}. 
Adaptation of our heuristic with the Siegel-Walfisz Theorem yields local bias factors 

Zf(p) - 1 

(55) b f , ai >0a (p) = <j P-*/(p) if € ^ . 

k ^/(p) if p e? 2 

Note that we still have b x k aiCl2 (p) = 1 for all p 6 "P, fc G N. 

Indeed, one can establish the following generalization of Theorem [TJ 

Theorem 3. // ai,a,2 £ Z, ai > 0, 02 7^ 0, f(x) = CkX k + ■ ■ ■ + c\x + cq G Z[x], > 0, and k > 1, then 

(ah TT fi 1 */(p)~M TT p~*/(p) TT - M ^h^l k N(k+1)/k 
/ 56 ^ pevo v wis peVl f p€V2 c k 1U & 

/JV( fc+1 V*l glogJV" 

+ 5 

V log 3 TV 

where the implied constant depends on f,a±, and a 2 . 

The vast majority of the argument for Theorem [3] is immediately analogous to the proof of Theorem [TJ so 
rather than starting from scratch, we proceed by highlighting the key differences. Much like (p~6|) . we start 
by giving a weighted, transform side formulation of our count, defining 

(57) Rf,ai,a 2 = / A N (aia)A N (a2a)SM(a)da, 

Jo 

where now 

max{ai, a\ + a 2 }A^ 1 



(58) M = 



Ck 



The procedure for absorbing the minor arcs into the error term goes through just as before, but some 
subtlety arises in the major arcs. For example, if a is in some major arc M / g , then ([2"B")) still gives an 

estimate for Ajv(aia), just potentially with a different denominator, based on the factors shared by ai and 
the original denominator q. More officially, we decompose each squarefree number based on our previous 
partition of the primes, i.e. we define 

(59) Nj = {pip 2 •■•Pf-Pj € Vi distinct} 

for i = 0, 1,2, where each set includes 1, so every squarefree number can be written uniquely as qm\m2 
with q G No, mi G Ni, and m 2 G N 2 . (Recall that due to the presence of the Mobius function, only 
major arcs with squarefree denominators contribute to the main term.) Now, one can check that by (|26[) . if 
a = ajqm\m,2 + (3, qm\m2 < log A, {a,qm\m2) = 1 and |/3| < log A/A, then 

(60) A^(a ia )7W(a 2 a) = ^}f ) .^}j ni \ v N {a l p) VN {a 2 P) + 0(A 2 e - c v^) 

<p{qy<j)(mi) 

for some c = c(S, ax, a 2 ) > 0. 

In place of the integral observed in ([34]) , we now have 

(61) / ( [ e- 2wia ^dy] ( ( N e- 2via ^dz ] ( [" e 2mc " x ^ dx] dp, 

J\/3\<\og B N/N \Jo J \Jo J \Jo J 

which with similar tricks of changing variables and applying Parseval's Identity can be seen to equal 
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(63) 



In place of the singular series ©(/), we now have 



&(f,a 1 ,a 2 )= ^2 



<t>{q) 2 4>(mi) qmim 2 



E 



'a/qmim2 : 



m 2 eN 2 



0<a<qmirri2 
(a 7 qmini2) — l 



which is still an absolutely convergent sum of a multiplicative function that is zero for non-squarefree numbers, 
hence we still have the Euler product formula 

6(/,ai, a 2 ) = Y[(l + F(p)), 
per 

where now 

1 



(64) 



F(p) = I 



PiP - I) 2 to 
1 

It,** 



a/p 



/ p if p e Pi • 
if p e -P 2 



Recalling our previous calculation that 



the result follows. 



p-i 

a=0 



□ 



Remark. Generalization of this result to linear combinations of more than two primes is definitely possible, 
but not particularly useful, as the circle method machinery is already capable of tackling the strictly harder 
task of counting solutions to a\p\ + • • • + aipi = N for fixed N and £ > 3. A useful pursuit, however, may 
be to extend the result to count solutions of the form a\p\ x + a,2p\ 2 = f(d), where fci and k 2 are arbitrary 
fixed natural numbers. 
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